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Abstract 

The Lindhard-Winther (LW) equipartition sum rule shows that within the 
hnear response theory, the stopping power of an energetic point-charge projectile 
in a degenerate electron gas medium, receives equal contributions from single- 
particle and collective excitations in the medium. In this paper we show that 
the LW sum rule does not necessarily hold for an extended projectile ion and 
for ion-clusters moving in a fully degenerate electron gas. We have derived a 
generalized equipartition sum rule and some related sum rules for this type of 
projectiles. We also present numerical plots for He"*" ion and He"*" ion-clusters. 
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1 Introduction 

The stopping power (SP), which is a measure of energy loss of energetic charged 
particles in a target medium, is of continuing theoretical and experimental interest 
in diverse areas such as interaction of charged particles with solids (see [1-4] 
for reviews) and beam- heating of fusion plasma For high- velocity projectile 
particles or clusters, the energy loss may be mainly due to collective and single- 
particle excitations in the target medium. It is of fundamental and practical 
interest to study the extent to which the collective and single-particle excitations 
each contribute to SP. This is the objective of our work reported here. 

The energy loss of a high velocity projectile is formulated, originally due 
to Lindhard 0, on the justifiable assumption of a weak coupling between the 
particle and a target medium which is modelled by a degenerate electron gas 
(DEG), through the linear response function of the DEG. The corresponding di- 
electric function e{k,uj) contains contributions from both collective excitations 
(plasmons) and single-particle excitations. For a single point-ion projectile Lind- 
hard and Winther (LW) investigated the respective contributions of these two 
excitations and found a sum rule which states that both these excitations con- 
tribute equally to SP. To our knowledge this type of sum rule has not yet been 
studied for an extended-charge projectile or a cluster. In this paper we formulate 
a generalized stopping power sum rule, again in the linear response approach, 
and present mostly analytical results for a He"'" ion and for a diproton cluster as 
projectiles. We compare and contrast our results with those of LW. 

In linear response theory, the stopping power (SP) which is the energy loss per 
unit length for an external projectile with a spatial charge distribution pext(i", t) = 
Qext (r ~ Vt) moving with velocity V in a homogeneous isotropic medium char- 
acterized by the dielectric function e{k,u!), is given by (see, e.g., 0]) 



where G(k) is the Fourier transform of the stationary charge (5ext(i")- 

Eq. (1) is applicable to any external charge distribution. In Sec. II we 
discuss a dicluster of two identical He"*" ions separated by a variable distance R. 
The spatial distribution of bound electrons in the ions, p(r), is assumed to be 
spherically symmetric. We use a Is-type wave function of the form '?/'is(r) = 
{Z^ /TTaQ)^^"^ e~^^^°'° , to describe the bound electron on each He"*" ion, with = 
0.529 A as the Bohr radius. Ze is the charge on each of the point-like nuclei 
{Z = 2). It may be remarked that we are considering an unscreened Is electron. 
For the projectile systems under study we may write G(k) as 




(1) 
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parameter. 

For a dicluster of He"*" ions the SP of a dicluster can be written as 



S = 2S'ind(A)+2S'eorr(A,i?,^) (3) 

— ^ / udu / Z^(z)lm— -\1 + COS (Az COS Jo (Bz sin {})]zdz^ 

7rA2 Jo Jq e{z,u) 



where A = {2u/ X) kpR, and B = 2kFR^yi — /\^. ^ind (A) and S'corr (A, -R, ^9) 
stand for individual and correlated SP, respectively. Jq{x) is the Bessel function 
of first kind and zero order and d is the angle between the interionic separation 
vector R and the velocity vector V; A = V/vp, = l/nkpaQ = (4/97r^)"'^^^ r^, 

1 /3 

Tg = (3/47r?7,oaQ) . is the electron gas density, vp and kp are the Fermi 
velocity and wave number of the target electrons respectively. In our calculations 
X (or Tg) serves as a measure of electron density. Here, as in Refs. [2-4,6,7], 
we have introduced the following notations z = k/2kp, u = uj/kvp. In these 
variables Z{z) = Z — p{z) and the Fourier transform of the spatial distribution 
p(r) is expressed as 
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where a = nx^Z . 

In Eq. (3) the term for correlated stopping power S'corr vanishes for large R 
[R — >■ oo) and SP is the sum of individual stopping powers for the separate ions. 
For i? — >• the two ions coalesce into a single entity. Then S'corr = S'ind and SP is 
that for a total charge 2e(Z — 1). 

We will consider the interaction process of a dicluster in a fully degenerate 
(T = 0) electron gas. For this purpose we use the exact random-phase approxi- 
mation (RPA) dielectric response function due to Lindhard 0, 



n=l 



2, ,2 X N r 2 

P + 2k ■ k„ - {2m 1%) (cu + zO+) 
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P - 2k ■ k„ + (2m/n) (tu + iO+) 



Here, k„ is the wave vector of the electron in the rath state. The distribution 
function / (/c„) is an even function of k„, and normalized so that iV = / 
is the total number of electrons. In the case of a fully degenerate free electron 
gas with Fermi energy Ep = h'^kp/2m the distribution function is / = 1 for 
kn < kp, and / (A;„) = for kn > kp. 

The summation in Eq. (5) can be analytically performed, leading to the 
characteristic logarithmic expression in e{z,u) first obtained by Lindhard H. 
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However, in our further consideration it is convenient to use the form (5) for the 
dielectric function. This allows a wider investigation of the analytical properties 
of e{z,u) in the complex z plane with applications toward deriving some useful 
sum rules. 



2 Equipartition sum rule 

In this section we shall discuss some important properties of the SP integrals, 
which we call stopping power summation rules (SPSR). In the literature two 
such sum rules have been widely considered— the familiar Bethe sum rule and 
the Lindhard-Winther (LW) equipartition rule (see 0). While the Bethe sum 
rule concerns the integral of l/e{z,u) over uj (or u) for fixed k (or z), the LW 
equipartition rule concerns the integral of l/e{z,u) over z for fixed u. The latter 
summation rule states that an integral of the form 

-1 

Im— -zdz = ^p(u) + ^spiu) (6) 



receives equal contributions, '^p{u) and '^sp{u), respectively, from the plasma 
resonance (plasmons), < 2; < m — 1, and from the region of close collisions 
(single-particle excitations), u — l<z<u + l. 

The LW equipartition rule was originally formulated for a single point-like 
charged projectile. Here we shall examine and generalize this equipartition rule 
for extended projectiles and their diclusters. As an example we shall consider a 
dicluster of two He"*" ions separated by a variable distance R. In order to deal 
with extended projectiles, as a generalization of Eq. (6), we need to consider an 
integral proportional to that in Eq. (3), 

f°° —1 

/ Z^{z)lm— -[l + cosiAzcosd) Jo (Bz sin d)] zdz, (7) 

Jo ^(^5 ^) 

which we decompose as '^p{a,u) + 53sp(a,'u), where '^p{a,u) and 53sp(a,'u) are 
the contributions for plasmon and single-particle excitations. 

In Eq.(7) the respective contributions for plasmon and single-particle excita- 
tions can be written as the sum of individual (first term) and correlated (second 
term) stopping terms 



%{a, u) = Qf\a, u) + %°"{a, u), Qsp{a, u) = %f{a, u) + 53™"(a, u). (8) 

Here again the terms for correlated stopping contribution, and Q'gp", for 

i? — >• tend to Q'"'^ and Q^^p respectively. 

Our objective is to show that the LW equipartition rule, '^p = Qgp, is not 
necessarily satisfied for extended projectiles and their clusters. We may therefore 
introduce the function 



4 



A(q;, u) 



(9) 



which represents the difference between single-particle and plasmon contributions 
to the integral given by Eq. (7). In order to calculate the function A{a,u), it is 
imperative to consider the integral in the complex z plane and to find the poles 
of l/e{z,u), i.e. the zeros of e{z,u), for fixed u. It is seen from Eq. (5) that 
there must be AN + 2 zeros of e{z,u) for fixed u. For a large real value of u, 
above the value Uq which corresponds to the minimum in the plasma resonance, 
all 4N + 2 zeros lie on the real z axis. Two zeros occur at 2; = ±Zr{u), and are 
determined by the intersection with the plasma resonance curve. 2N zeros are 
grouped together in the intervals u — 1 < \z\ < u + 1, and are responsible for the 
single-particle contribution, '^gp, to the integral in Eq. (7). The remaining 2A'" 
zeros are also grouped together near the point z — 0, but they he outside of the 
interval u — 1 < \z\ < u + 1 which is responsible for single-particle excitations. 
It can be seen directly from Eq. (5) that at 2; ^ 6{z,u) — ^ 00. Therefore the 
contribution of these latter 2N zeros of e{z, u) to the integral in Eq. (7) vanishes. 

It should be noted that in the RPA dielectric function (5), the variable u has 
an infinitesimal positive imaginary part which is introduced for causality. It is 
seen from Eq. (5) that a small positive imaginary part iO"*" being added to u, 
is equivalent to a displacement of the zeros of e{z^u) in the complex z plane in 
such a way that the zeros z = ±Zr{u) at plasma resonance lie below the real z 
axis while the zeros in the region u — 1 < \z\ < u + 1, lie above this axis. Prom 
Eq. (5) it also follows that the slope de{z, u)/dz is positive at the plasmon curve, 
z = Zr{u) but negative at the other zero curves, z ~ Zj{u) in the z plane. 

With these observations we shall now establish a generalized SP equipartition 
sum rule for a He"*" ion di cluster projectile, and shall present results for two values 
of the orientation angle i? made by the projectile velocity vector with the inter- 
He+ ion separation vector. The result for a single He+ ion will subsequently be 
obtained from the dicluster sum rule. 

2.1 He''" ion cluster with 'd = 

To derive an expression for a non-zero A (a, u) and to generalize the LW summa- 
tion rule we consider the following contour integral 



where the contours Ci (for Q^{a,u)) and Di (for Q_{a,u)) are shown in Pig. 1. 
These two contours contain the real z axis (0, -l-oo), upper (for Ci) or lower (for 
Di) quarter circles, the imaginary z axis (±ioo,0) and infinitesimal semicircles 
C2 or D2. The full structure of the integral in Eq. (10), for the cases of He+ ion 
and a dicluster, is more involved than for the case of a point-like ion, and it also 




[1 -|- exp {±iAz)] zdz, 



(10) 
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contains an exponential function. An analytical evaluation of this integral thus 
leads us to consider a contour different from the one used in the LW paper (see 
1^, for details). 

For large values of the dielectric function must behave as 



e{z, u) 



1 + — 



according to Eq. (5). Therefore both integrals Q±{a,u) vanish within the upper 
and lower quadrants, respectively and from Eq. (10) we find 



Q±{a,u) 



J±ioo JC2,D2 



zHz) 



e[z,u) 



[1 + exp (±i^2;)] zdz. 

(12) 

On the other hand both the integrands in Q±{a, u) are analytical functions in- 
side contours Ci and Di containing single-particle N poles Zj+iQ^ , j = 1,2, ...,N, 
(contour Ci) or single plasmon pole, Zr{u) —iO~^, (contour Di). According to the 
theorem of residues, for these functions we have 



Q+{a,u) 



-2m ^ ^^^\ [1 + exp [iAzj)] 



Q^{a,u) = 2ni ^^ ^^''^ 

Q^ByZr, U) 



[I + exp (iAzr)] 



(13) 
(14) 



As has been mentioned earlier, de{zr,u)/dz > for plasmons and de{zj,u)/dz < 
for single-particle excitations. 

Now let us take the imaginary part of both sides of Eqs. (12)-(14). We find 



I±{a,u) 



Im / Z\z) 

'0 



Q'p(a, u) 



where 

P±{a,u) = : 
In Eq. (15) we have introduced 



[1 + exp {±iAz)] zdz (15) 



P±{a,u), 



[ Z\z) 




'C2,D2 


e[z, u) 



[1 + exY>{±iAz)]zdz. (16) 



q\ ^ ^[ [1 + COS {Azj) 



Z'^^z) Im- 



[1 + COS {^Az)\ zdz. 



(17) 
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u-1 



e{z,u} 



[1 + cos (Az)] zdz, 



for single-particle (53 ^^(q;, u)) and plasmon (3p(a;, u)) contributions to the integral 
(7), respectively. To prove that the sums in Eqs. (17) and (18) are actually equal 
to the integral forms of Q'sp(Q;, u) and Qp{a, u) we use the known expression (see, 
e.g., i) 



Im 



e{z, u) 



TT 



Im e(2,u)— >0+ 



5 {z — Zj) 



te[zj,u) 



(19) 



where Zj are the zeros of s{z, u 



Since all single-particle poles lie in the interval 
u — l<z<u + l, the integration in Eq. (17) over z results in the summation 
form for Qsp{<^,u). However, for plasmon contribution the plasmon pole Zr{u) 
lies in the interval < z < u — 1, only for sufficiently high u, u > uq{x) 
The threshold value 'Uo(x) depends on the electron gas density. For instance, for 
metallic densities ~ 2 (x ~ 0.5) for threshold value of u we have Uq ~ 1.4. 
When u < uo{x) the plasmon contribution term vanishes, Q'p = 0. 

Next, we note that the imaginary part of the second terms in Eq. (12) has 
been omitted, because the function 



Z'(±iz) 



1 - 



1 



e{±iz, u) 



(20) 



is real. Consequently 



Im 
Im 



zHz) 



±ioo 
oo 



1 - 



e{z, u) 



Z'(±iz) 



e(±iz, u 



[1 + exp {±iAz)] zdz (21) 
1 + exp (—^2;)] zdz = 0. 



Let us now consider the integrals P±{a,u). From Fig. 1 and Eq. (16) and 
after evaluating the residues we find 



P+{a,u) 



P-{a, u) 
3 



(22) 



1 

16 



d a d'^ 
+ 



a 



Q3 



da 16 48 da^ 



where 
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= [1 + exp {-Aa) 



1 



e[a, u) 



(23) 



lij = e\ia, u) — 1 -\- - — r < —2a — arctan h arctan (24j 

Aa"^ y a a 

t 1 i\ ( Oi a \ ^ (u + l)^ + a^} 

+ (u — a ) [ arctan arctan +ua In — > . 

^ ^ \ u+1 u — lj [u — iy + a-^) 

Therefore from Eqs. (15) and (22) we finally obtain 

IS.f\a,u) = P^{a,u) + ]-^+{a,u)-I-{a,u)] (25) 



2 



j'OO 

= P+(a,M) + P / Z'^{z)Re 
Jo 



e{z, u) 



sin (Az) zdz. 



Note that in Eq. (25) a,t z — Zr{u), where < Zr{u) < u — 1, both the real and 
imaginary parts of the dielectric function vanish i.e. e{zr, u) = 0. Therefore one 
needs to consider the Cauchy principal value (denoted by P) of the integral. 

Eq. (25) is the generalized SPSR for a lle+ ion cluster with '& = 0. From the 
general expression (25) we shall derive below some particular SPSR for individual 
He+ ion and diproton cluster with — 0. 

2.1.1 Individual He"*" ion 

We can derive the SPSR for an individual He+ ion directly from Eq. (25) if 
we consider the limit oi A ^ oo. In this limit the exponential function in Eq. 
(23) eventually vanishes and the function P^{a,u) is defined by the first term 
in Eq. (23). The second term in Eq. (25) in the limit of ^ — ^ oo must behave 
as cos [a (m) ^] as can be seen from Eqs. (17) and (18), where a{u) is some 
unspecified function of u. Therefore the second term in Eq. (25) oscillates with 
an increasing A or interionic distance R; the full integral of this term over u is 
damped although not necessarily vanishing as .4. — oo. However, when we include 
a small damping in the electron gas, which is expected for any real medium, the 
second term in Eq. (25) vanishes as ^ — > oo. Thus for an individual He^ ion we 
find: 

A^''\a,u)^Poo{a,u), (26) 

where P^{a,u) is the function P^{a,u) at .4. oo. 

For individual protons {a oo) the right-hand side of Eq. (26) behaves as 

P^(a,«)^|^(Z-l)^0. (27) 
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Consequently, in this limit we recover the known Lindhard-Winther equipartition 
rule (ER), = %. 

2.1.2 Diproton cluster with ^ = 

For a diproton cluster (a oo) the function P+{a,u) vanishes as in Eq. (27). 
Therefore for a cluster of point-like particles we find 



For individual protons (i? — * oo or ^ — > oo) the RHS of Eq. (28) vanishes 
due to a small damping in the electron gas and again we recover the Lindhard- 
Winther ER. 

2.2 He+ ion cluster with = 7r/2 

In order to derive an analytical expression for A{a,u) for a He"*" ion dicluster 
with 19 = 7r/2, we use the same integration contours Ci and Di of Fig. 1, and 
consider the following integrals 



where Hq (z) and Hq (z) are the Hankel cylindrical functions of the first and 
second kind, respectively and of zero order. We may recall that the Hankel 
functions are analytic inside and on the contours Ci (for Q^{a,u)) and Di (for 
Q-{a,u)) except at the point z = where they have a logarithmic singularity. 
Moreover, the functions H^^^ (z) and H^^^ (z) vanish on the upper and lower 
quadrants, respectively. 

Using the theorem of residues, the functions Q±{a,u) are evaluated as 




(28) 




1 + i^z) zdz, 



(29) 




(30) 



Then for i9 = tt/2, we obtain from Eqs.(29)-(31) 




1 + (Bzr)] 



(31) 




(32) 



$^p(q;, u) 



7lZr2'^{Zr) 



[1 + Jo (BZr)] , 



(33) 



^ a 
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which are different from the corresponding quantities for ^ = given in Eqs.(17) 
and (18). 

Using a similar procedure of calculation as in the previous section we finally 
find 



AW(a,«) 



Poo (a, u) + F 



TT 



POO 


1 


/ Z\z) Re 




Jo 


e{z, u) 



Yo{Bz)zdz (34) 



POO 




/ 2\z) 




Jo 


e[z, u) 



Kq {Bz) zdz, 



where Z{z) = Z{iz), Yq {z) is the Bessel function of the second kind and zero 
order, Kq {z) is the modified Bessel function of the second kind and zero order. 

After some algebraic manipulation, the last term in Eq. (34) can be shown 
to read: 



TT 

+ 

TT 



e{z, u) 
1 \ d 



16 / da 



Kq {Bz) zdz 



+ 



a (9^ a^ 


p r Kq {Bz) 


1 


16(9a2 48 9a3_ 


Jq a^ - z^ 


e{z, u) 



(35) 
zdz. 



Eq. (34) is the generalized equipartition sum rule for a He^ dicluster with 
= 7r/2. In the limit a ^ cxd, we obtain the corresponding sum rule for a 



diproton cluster with t? = 7r/2, 
dSP\a,u) = 



POO 




I Re 




Jo 


e{z, u) 



2 

TT 



POO 


1 






Jo 


e{z, u) 



Yq {Bz) zdz 
Kq {Bz) zdz. 



(36) 



The first and second terms in Eq. (36) vanish as B^^^"^ and B~^ respectively 
at large interionic distances. This leads again to the LW equipartition rule for 
point-like projectiles. Thus, in contrast to the case of an aligned dicluster with 
"i? = (see Eq. (28)), for d = 1^/2 there is no need to introduce an infinitesimal 
damping for plasmons and single-particle excitations to get the correct limit at 
B ^ oo. This is because a fast ion moving in an electron gas excites the wake 
field which has different structures along and across of the direction of motion 
13] . These waves are damped strongly across the direction of motion P|. 



2.3 Some simple SPSR 

In this subsection we will consider some useful summation rules. To derive the 
first of them we consider the following contour integral 
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Q{a,u) 



[ Z\z) 




z In zdz — — 27ri 




e[z^ u) 





■ In Zj, (37) 



where the contour Ci is shown in Fig. 1. 

Now, in contrast to the previous sections we take the real part of both sides 
of Eq. (37). Then using similar calculational techniques as in Sec. 2.1 we find 



poo 

Jo 



Re 



e(z,u) 



TT 



z In zdz 

1-^ 



(38) 



z In zdz. 



In a similar way we can obtain some other summation rules 



/■oo 

' / Z''{z) Re 
Jo 



s{z, u) _ 



zdz 



^0 



Z\z) 



U) 



zdz, (39) 



and 



POO 


1 


/ Re 


e{z,u)_ 


Jo 



z In zdz 



POO 


1 
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e{z,u)_ 



z In zdz 



(40) 



Re 



1 



e[z,u] 



zdz — — I 1 — zdz. (41) 

Jo V £{z,u)\ 

The singularities in the RHS of Eqs. (38) and (39) can be understood in a similar 
manner as in Eqs. (34) and (35). 

In order to illustrate our analytical results, Eqs.(25)-(28) and (34)-(36), in Fig. 
2 we present A(q;,-u) as a function of the parameter u = uj/kvp {uq ^ -u ^ A, 
where Uq is the threshold value for plasmon excitation) for an individual He+ ion 
(solid line). The lines with and without circles correspond to He+ ion (R = 3 A) 
and diproton (R — 1 A) clusters, respectively, = (dashed hues) and ■& — 7r/2 
(dotted lines). The numerical results are for fast projectiles X — V/vp — S and for 
density parameter — 2.07 appropriate to the valence electrons in Al. A{a,u) 
as a function of u has interesting features. A{a,u) remains positive {'^sp > '^p) 
for He"*" ion and for He"*" ion cluster with = (at m > 2) and = 7r/2 and 
remains negative C^gp < '^p) for diproton cluster with i? = and ■& — The 
curve for He+ ion cluster with i? = crosses the zero-axis for w ~ 2. From Fig. 
2 the oscillatory nature of the function A(q;,m) can be seen for both diproton 
and He+ ion clusters with -i? = 0. Note that for a diproton cluster with d = tt /2 
at M = A, when the excited wave moves with the phase velocity uj/k = V, the 
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parameter B vanishes, B — 0. Therefore, from Eqs. (36), (40) and (41) at B — > 
follows A = 0, i.e. the LW equipartition strongly holds. While for He"*" ion cluster 
with i} = 71 /2 from Eqs. (34), (38) and (39) we find that the function A at m = A 
is two times greater than for an individual He^ ion projectile. In general for 
the high velocity domain uj/k ^ vp {ox u ^ 1) all the curves decrease and an 
approximate LW equipartition rule, ~ '^Sp holds asymptotically. But it is also 
clear that in this high velocity limit the energy losses due to single-particle {^sp) 
and plasmon ($5p) excitations decrease as well. 

We may conclude with the remark that our analytical expressions are well 
supported by numerical results. 
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Figure Captions 



Fig. 1. Illustration of contours Ci and Di, in complex z plane. Isolated point 
P below real z axis indicates plasmon pole. Group of crosses above real z axis 
indicates poles in single-particle excitations. 

Fig. 2. The function A(a,u) vs parameter u for A = 8. Individual Hc+ ion 
(solid line), diproton cluster with R = 1 A, i9 = (dashed line without circles) 
and d = Ti/2 (dotted line without circles), He^ ion cluster with i? = 3 A, -i? = 
(dashed line with circles) and 'd = ti/2 (dotted line with circles). The density 
parameter is Vg — 2.07 (Al target). 
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